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We consider shock probes in a one-dimensional driven diffusive medium witli near- 
est neigfibor Ising interaction (KLS model). Earlier studies based on an approximate 
mapping of the present system to an effective zero-range process concluded that the 
exponents characterising the decays of several static and dynamical correlation func- 
tions of the probes depend continuously on the strength of the Ising interaction. On 
the contrary, our numerical simulations indicate that over a substantial range of the 
interaction strength, these exponents remain constant and their values are the same 
as in the case of no interaction (when the medium executes an ASEP). We demon- 
strate this by numerical studies of several dynamical correlation functions for two 
probes and also for a macroscopic number of probes. Our results are consistent with 
the expectation that the short-ranged correlations induced by the Ising interaction 
should not affect the large time and large distance properties of the system, implying 
that scaling forms remain the same as in the medium with no interactions present. 



I. INTRODUCTION 



Useful information about a complex system is often obtained by introducing probe parti- 
cles into it. After the probe particles have come to a steady state with the system, their static 
and dynamic behavior often reflect important characteristics of the system. For instance, 
by monitoring the motion of probe particles, one can understand visco-elastic properties of 
a cel. II the s„.-ge. tran.ifo,. in a polymer solution or eorrelations present .n bacterial 
motion p|. In certain cases, for example, in active micro- rheology, the probe particles are 
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subjected to external force fields. In the forced dynamics of a magnetic bead in a dense 
colloidal suspension, has been used to study the colloidal glass transition. 

In this paper we will consider one such example of nonequilibrium (driven) probe particles 
which are introduced in a nonequilibrium medium, to study how the static and dynamic 
properties of the probe particles are influenced by the surrounding medium, and also how 
the medium gets affected by the presence of the probe particles. 

We study a particular simple one-dimensional lattice gas model first introduced in [sl to 
describe the motion of probe particles in a current-carrying medium. The probe particles are 
taken to exchange with particles and holes of the medium with equal rates but in opposite 
directions. Because of these dynamical rules, the probe particles tend to migrate towards 
the region of strong density variations (or shocks) which may be present in the system. 
Studying the dynamics of these shock-tracking probe particles, one can therefore infer the 
motion of density fluctuations in the medium. 

In an earlier study [6|], we have discussed the dynamical properties of these probe parti- 
cles in a nonequilibrium current-carrying medium in which there is no interaction between 
medium particles except hard-core exclusion. In this case, the medium was described by an 
asymmetric simple exclusion process (ASEP) which is the simplest lattice model of driven 
diffusive systems The shock-tracking probe particles then reduce to second class parti- 
cles [Sj]. Derrida et al. have found the exact stationary measure of the system [9]. Their 
studies on static properties of the system show that when the number of second class par- 
ticles is finite, they form a bound state and the steady state distribution function of the 
separation r between a pair decays as r^^ with A = 3/2. A macroscopic number of second 
class particles gives rise to a correlation length which diverges (proportional to the square 
of the interprobe separation), as the probe concentration goes to zero. 

We studied the dynamical properties of this system in presence of a macroscopic number 
of probe particles, and found that the dynamics is governed by a time-scale which marks 
the crossover from single-probe behavior to many-probe behavior [g]. This time-scale shows 
a strong divergence (proportional to the cube of the interprobe separation) in the limit of 
vanishingly small density of the probe particles. This diverging time-scale is related to the 
diverging correlation length present in the system y], and enters the scaling descriptions of 
various dynamical correlation functions of the probe particles jo]. 

In the present paper, we present a detailed study of shock-tracking probes in a driven 
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system in which there is a short-ranged Ising interaction between the particles of the medium. 
In the absence of any probes, such a medium can be described by the 1-d Katz-Lebowitz- 
Spohn (KLS) model, whose steady state has an Ising measure 10, 11|. In [sl Kafri et al. 
have reported that in presence of a macroscopic number of probe particles, the system 
shows an interesting phase transition as the strength of the Ising interaction is varied. 
Beyond a critical value of the interaction strength and for sufficiently high density of the 
medium, a macroscopic domain consisting of particles and holes (no probes) is formed. A 
characterisation of this phase transition was attempted using an approximate mapping to 
the zero-range process (ZRP) where the probes are mapped onto ZRP sites and the particle- 
hole domain preceding a probe is mapped onto the occupancy of that site The current 
out of a particle-hole domain then becomes the hopping rate out of a site in the ZRP. 

A prediction of this approximate mapping is that exponents characterizing the decays of 
several static and dynamic quantities should depend continuously on the Ising interaction 
strength e. However, our numerical studies of these quantities seem to indicate that the 
exponents are e-independent over a substantial range of e. This paper is concerned with a 
study of the differences between our results and those based on the ZRP picture. 

A possible simple rationalization of our results is that the Ising interactions would be 
expected to give rise to a finite correlation length Rising, whose value may be renormalized in 
the presence of probes, but is still expected to be finite. Then, on length scales r ^ Rising, 
the system should behave essentially like the non-interacting (e = 0) system. Our results are 
indeed consistent with such a scenario, as we find e-independent behaviour asymptotically 
(for large r and t), even though there is sometimes an e dependence for smaller r and t. 

Below we describe in brief the quantities we studied and our results. 

(1) Distribution function of the size of the particle-hole domains: The mapping to the 
ZRP predicts that in the disordered phase this distribution function should be an ex- 
ponential times a power law, with a power which is a continuous function of the Ising 
interaction strength. However, we observe that the domain size distribution shows 
a power law exponent which does not vary with e but remains constant at its value 
for e = 0. To understand this discrepancy, we are led to check the assumptions that 
have been made in the approximate KLS-ZRP mapping. We find that the assump- 
tion of statistical independence of the domains remains valid, and further verify that 
accounting for the finite size correction to domain currents is not the reason behind 
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the discrepancy. However, as we discuss in section 3, the movement of probes in a 
KLS medium is non-Markovian and the ZRP mapping does not capture this aspect of 
probe dynamics. 



(2) Dynamics of two probe particles: This was studied in |12l.ll3l| where it was reported that 
a bound state forms between the probe pair such that the distribution of separation 
decays as a power law with an exponent 6(e) that varies continuously with the strength 
of the Ising interaction e. Starting from a configuration where the two probes were 
nearest neighbors, a scaling form was proposed to describe the temporal evolution of 
their separation. The authors have tried to verify this scaling form by measuring the 
cumulative distribution and the mean value of the separation between the probes as 
a function of time. They reported that in conformity with their scaling hypothesis, 
the time-dependent cumulative distribution function P{r,t) for different values of t 
undergoes a scaling collapse when rescaled by and plotted against rt~^^^, 
where z is the dynamical exponent that takes the value 3/2. The average separation 
between the two probes is reported to grow with time as a power law with an exponent 
[2 — b{e)]/z, which is consistent with their scaling form. 

On the contrary, we find that although for an initial time range, the average distance 
does show e-dependent growth, for larger times, it crosses over to another growth 
regime where the exponent takes a value which is close to the one expected for e = 0, 
i.e. with no Ising interaction. Our numerical results also show that for larger times, 
the scaling collapse of P(r, t) fails. We have verified that the scaling collapse can be 
retrieved by rescaling with t^^^ (instead of jf:[^(^)^i]/^)j as in the case of e = 0. 

(3) Dynamical properties with a macroscopic number of probes: The dynamical correlation 
functions in this case are found to follow the same scaling description as with e = 0, 
with a crossover time-scale which separates a single-probe regime at short times from 
a long-time regime characterized by collective behaviour of the probes. Moreover, 
the crossover time-scale shows a similar divergence in the limit of vanishingly low 
concentration of the probe particles. In other words, our studies indicate that even in 
the presence of a nearest neighbor Ising interaction in the medium, the large time and 
large distance properties of the system do not change. 

In the following section, we describe the lattice model on which we have performed Monte 
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Carlo simulation and briefly summarise our earlier results for the non-interacting (e = 0) 
medium. In section 3 we discuss the static properties of this model where we recall the 
approximate mapping to the zero-range process (ZRP) introduced in js] and discuss the 
validity of various assumptions that went into this mapping. In section 4 we discuss the 
dynamical properties of the system in presence of a flnite number of probes and also for a 
flnite density of the probes. 

II. MODEL AND EARLIER RESULTS (e = 0) 

The model is deflned on a one dimensional periodic lattice each site of which may either 
be empty or may contain a particle of the medium or a probe. We use the symbol '+' to 
denote a particle, '— ' to denote a hole and '0' to denote a probe. The exchange rules are as 
follows. 

+0^0+ (1) 
0- -0 

Here AV^ is the change in the nearest neighbour Ising interaction potential 



V = -jY,s,s,+, (2) 
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where Si = 0, ±1, according as the site i contains a probe, a particle or a hole, respectively. 
Throughout we consider equal densities of particles and holes in the medium, i.e. po = 1 — 2p 
where p and po denote densities of particles and probes, respectively. The coupling parameter 
e may vary in the range [—1, 1]. In this paper, we will only consider e > 0. 

In the absence of any probes, the system reduces to a 1-d KLS model with an Ising 
measure in the steady state. This gives rise to an e-dependent correlation length Rising in 
the system. For e < 1, this correlation length remains flnite and hence the large distance 
properties of the system can be expected to remain unaffected by the interaction. 

When probes are present, as seen from the last two exchange rules in Eq. [H a probe 
exchanges with particles and holes of the medium in opposite directions. This implies that a 
probe would tend to be located in a position where there is an excess of holes to its left and 
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an excess of particles to its right. In other words, there would be a strong density variation 
or 'shock' around a probe. This is the reason we call them 'shock-tracking probes' (STPs). 

In the absence of any interaction, one has e = and in this case, a particle in the medium 
executes a totally asymmetric exclusion process (TASEP) with an effective hole density 
(1 — p) and it exchanges with a hole and a probe in the same way; similarly, a hole in the 
medium also executes a TASEP [in the opposite direction and with an effective particle 
density (p + po)] and it exchanges with a particle and a probe in the same way. In other 
words, for e = a probe behaves like a particle for an adjacent hole, and like a hole for an 
adjacent particle, ^uch probes are known as 'second class particles' [8|]. 

Derrida et al. [3] have found the exact steady state measure of this system of second class 
particles in an ASEP by using the matrix method. In presence of more than one second 
class particle, the steady state factorises about any second class particle, which implies 
factorisation in terms of the one component system about the shock position. When there 
is a single second class particle present in the system, the mean density profile around it 
decays as a power law with an exponent 1/2. In the presence of two (or a finite number 
of) second class particles, the medium induces an attraction between them and they form 
a weakly bound state where the distance r between two successive second class particles 
follows a power law distribution P(r) ~ When the number of second class particles 

is macroscopic, the density profile at a distance r from any given probe takes the form 

P{r) ~ ^ exp (-r/ + p (3) 



where the correlation length ^ diverges in the low concentration limit of the probes 9|: 

e ^ 4p(l - p)/pI as Po ^ (4) 

We monitored several quantities to study the dynamical properties of systems with macro- 
scopic number of probes. We find a diverging time-scale which marks the crossover between 
single-probe behavior and many-probe behavior. In Section IV, we will discuss the behaviour 
of these quantities, when e is nonzero in the KLS model. 

The variance of the displacement of the tagged probes is defined as 

Co(t) = {(Ykit) - n(0) - - YkiO))?) (5) 



where yfc(t) is the position of the k-th probe at time t. Ferrari and Pontes l4[ had earlier 



calculated the asymptotic (t oo) behavior of Co{t) and shown that Co{t) ~ Dt with 
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diffusion constant D = [p(l — p) + (p + po)(l — P — Po)] / Po- For small times, in the limit 
of low concentration of the probe particles, one would expect each probe to behave as an 
individual non-interacting particle subject only to the fluctuations of the medium. The 
variance of the displacement of a single probe is found analytically to grow as 
In the limit of small but finite concentration of the probe particles Co{t) shows a single- 
particle (super-diffusive) behavior at small time and diffusive behavior at asymptotically 
large times. One would therefore expect a crossover between these two regimes that would 
occur at a time-scale r which is a function of po- The natural expectation would be r ~ 
where ^ is the correlation length as defined in Eq[3l Substituting the value of the dynamical 
exponent 2; = 3/2 and using Eq. Hlone obtains 



r ~ Po 3 (6) 



in the limit of small po. This leads us to propose the following scaling form for Co(t) 

C,^) - t"^F . (7) 

This form is valid in the scaling limit of large t and large crossover time-scale r (i.e. po — > 0). 
Here F{y) is a scaling function which approaches a constant as ?/ — > 0. For ?/ ^ 1, we must 
have F{y) ~ y^^^^, in order to reproduce Co(it) ~ Dt. We have verified the above scaling 
form by Monte Carlo simulation 6|. 

The same crossover time-scale r is found to be present in other dynamical correlation 
functions as well. To track the dissipation of the density pattern of the second class particles, 
we considered the quantity 



Bo{t) = {Yu{t) - n(0) - (n(t) - n(0))) (8) 

where the overhead bar denotes averaging over different evolution histories, starting from a 
fixed initial configuration drawn from the steady state ensemble (see IT], [iSj and also 6| 
for a discussion on why this special averaging process is useful in tracking dissipation). Our 
scaling analysis leads to the following scaling form: 



B,{t) ~ t^'^G (^-^ . 



(9) 



where r is the same crossover time-scale as in EqE] and G{y) is a scaling function which 
approaches a constant as ?/ ^ 0, while for 1/ ^ 1, one expects G{y) ~ y"^^^. Our numerical 
results are consistent with this scaling form 
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Finally consider the quantity 

A(t) = mt) - nm') (10) 

which measures how the separation between two successive probes fluctuates in time. Here, 
R{t) is the separation between the k-th and {k + l)-th pair at time t. Our studies show that 
A(t) has the following scaling form 

A(t)~tif(^) (11) 

where the scaling function H{y) approaches a constant as y — > and for y ^ I one must 
have H{y) ~ 

To summarise, for e = we flnd that several dynamical correlation functions of the 
probe particles are governed by a single crossover time-scale r which diverges as Pq ^ for 
low concentration of the probes. In the remaining portion of the paper, we will consider 
static and dynamical properties for e > and examine how different they are from the 
non-interacting case. 



III. STATIC PROPERTIES OF KLS MODEL WITH PROBES 

Kafri et al. reported that the KLS model with macroscopic number of probes shows 
phase separation transition for e > 0.8 as the density p is increased above a critical value 
Pc jsi. They concluded that in the phase separated state, a macroscopic domain, composed 
of particles and holes of the medium, coexists with another phase which consists of small 
domains of particles and holes, separated by the probes. They explained this phase transition 
by attempting to approximately map the system onto a zero-range process. 

To describe the mapping, we flrst deflne a domain as an uninterrupted sequence of parti- 
cles and holes, bounded by probes from both ends. The current J„ out of a domain of length 
n can then be determined by studying a KLS model in an open chain with boundary rates 
of injection and extraction equal to the rate at which the particles and holes of the domain 
would exchange with the probes at the domain boundaries. According to Eq. [1] this rate is 
unity. The current J„ can be calculated exactly for an open KLS chain and for large n it 
has the form 

Jn = JooU + ^] (12) 
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where the coefficient b has the following dependence on e 

The study of Kafri et al. indicates that b plays an important role in characterising the phase 
separation transition in the model. 

The present system is mapped onto a zero-range process (ZRP) as follows: the i-th probe 
is defined as the i-th site of ZRP and the length of the domain to the left of the i-th probe 
is taken to be the occupancy n{i) of the i-th site of ZRP. We illustrate this in figE The 
density in the ZRP is related to the KLS model density p as pzRP = 2p/(l — 2p). 



+ + -0 + - + -- 



FIG. 1: A typical configuration of the KLS model with probes and its corresponding configuration 
in ZRP. 



The hopping rate out of the z-th site in the ZRP is taken to be the domain current Jn(i) 
given in Eq. [121 For such a ZRP, the condition for condensation to take place is 6 > 2 
and pzRP larger than a certain critical density pc- In the condensed phase, the occupancy 
at a single site becomes macroscopically large, while the remaining sites have an average 
occupancy pc [IS!]. For pzRP < Pc, the number of particles present on a site follows the 
distribution function 

P{n) ~ ^ exp{-n/^zRp) (14) 



n'' 

where the correlation length ^zrp diverges as pzpp — > Pc- For pzpp = Pc, we have P{n) ~ 



1 



while for pzpp > Pc, a similar power law decay describes the distribution at all sites except 
for the single condensate site. 

The approximate ZRP correspondence implies that in the KLS chain with probes, for 
large enough p and for e > 0.8 (as follows from Eq. [13]), there should be a macroscopic 
domain present in the system which is composed of particles and holes (no probes). The 
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rest of the system should consist of small probe clusters, interrupted by the domains (of 
particles and hole) with size distribution given by Eq. [TH 

From numerical simulations for e < 0.8, it was found that a very large domain may exist 
for large p [5] . Our numerical simulations confirm this. In j^, it has been argued that 
this is not true phase separation, but rather a consequence of the fact that the correlation 
length in this case has a large (but finite) value. 

According to the above correspondence with the ZRP, it is expected that close to the 
critical point, the domain size distribution for n <^ ^zrp should follow a power law with ex- 
ponent 6(e) which should increase monotonically with e. However, our numerical simulations 
for various values of e and p [see fig |2] show that the power law exponent seems throughout 
to be close to 3/2 (which is the value of 6 at e = 0), independent of the value of e. 
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FIG. 2: Domain size distribution P{n) scaled by v?!'^ shows a flat stretch over a substantial range 
of n for different values of e and p. For comparison with the ZRP prediction, we have scaled P{n) 
by n^^^^ in the inset. In both these plots, the symbol + corresponds to e = 0.9, p = 0.464, symbol X 
corresponds to e = 0.8, p = 0.375 and symbol * corresponds to e = 0.6, p = 0.375. In the last case, 
€ is substantially smaller than the critical value and ^zrp is shorter. This explains the observed 
deviation from the power law behavior for large n. 
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This points to a contradiction with the correspondence with the ZRP, and leads us to 
examine the assumptions that go into the KLS-ZRP mapping. 

Independence of Domains: A crucial property of the ZRP is that the occupancies at 
the sites are uncorrelated. In our present model of the KLS chain with probes, this would 
imply that the domains between the probes should be independently distributed. We have 
verified this assumption by measuring the conditional probability P{n\n') that the size of 
a particular domain is of length n given that its neighboring domain is of length n'. We 
find that P{n\n') does not depend on n' and is same as P{n) consistent with neighboring 
domains being distributed independently. Our data is presented in fig [3l 
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FIG. 3: The conditional distribution of domain size P{n\n') as a function of n for n' = 4,8. For 
comparison P{n) is also shown. P[n\n') is seen to match with P{n) which shows the domains are 
independently distributed. We have used L = 2048, e = 0.6 and p = 0.375. 



Finite Size Correction to Domain Current Apart from the independence of domains, 
another requirement for the ZRP mapping to hold is that the current out of a domain of size 
n should be the same as the current in an isolated open KLS chain and is given by Eq. [121 
Evans et al. have given evidence for this by numerically measuring the actual current out 
of a domain and comparing with the exact calculation for an open chain KLS model 21|. 
Good agreement was found for large n. 
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To take into account the finite size corrections for moderate values of n, we simulated a 
ZRP where the hopping rate out of a site is read off directly from the actual J„ vs n data, 
obtained from numerical simulation of the KLS model with probes. The mass distribution 
for this ZRP is found to have the same form as in Eq. [14] with the exponent h given by Eq. 
[T3l as expected. We conclude that the finite size correction to J„ is not the reason for the 
discrepancy shown in fig [2j 

Non-Markovian Movement of the Probes: There is however, one aspect of the KLS model 
with probes that is not captured in the corresponding ZRP. Since a probe exchanges with 
the particles and holes of the medium in opposite directions, as shown in Eq. [H once a 
probe moves in one particular direction, it cannot move in the opposite direction at the 
very next time-step. For example, suppose a probe moves to the left by exchanging with a 
particle in the medium. Immediately after this exchange the probe has the particle as its 
right neighbor. Clearly, the probe cannot take a step to the right as long as that particle 
stays there. In other words, the probes have a finite memory which makes their movement 
non-Markovian. In terms of the ZRP this would mean that once a site has emitted a particle 
to its right neighbor, it has to wait for some time till it can receive a particle from its right 
neighbor. This waiting time should depend on the form of the density profile in a domain. 
Note that in this non-Markovian ZRP, apart from and 6(e), there are other parameters 
that are associated with the exact form of the waiting time. As a result, the phase-diagram 
becomes complicated and to specify the criterion of a phase transition a much more detailed 
analysis is required. This might shed some light on the observed discrepancy about domain 
size distribution. 



IV. DYNAMICS OF PROBES IN THE KLS MODEL 



Two Probes 



The properties of two STPs in the KLS chain were first studied in [12] by Levine et. ai, 
who argued that the time-evolution of the separation between the probe pair is governed by a 
Master equation. Their analysis indicates that the medium induces an attraction among the 
probe particles and they form a bound state. The steady state distribution of the distance 
between two probes takes the form P{r) ~ r^^ where 6 is a function of e given by Eq. [131 
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For e = one retrieves P(r) ~ r~^/^ as found in j^. 

Rakos et al. have shown that the random force between the probe pair is sensitive to the 



noise correlations present in the medium jl3|]. When the probe particles are embedded in a 
KLS ring, such that the random force that drives the probe particles is fully generated by 
the current fluctuations of the driven medium, the probes inherit the dynamical exponent 
of the medium, which is 3/2. On the other hand, if the random force has a part that is 
temporally uncorrelated, the resulting motion is described by a dynamical exponent z = 2. 

To study the dynamics of the system, the distance between the two probes was monitored, 
starting from the initial configuration in which the two probes were side by side. The 
approach to the steady state was modelled by the scaling ansatz 

P{r,t) ^ r-^f{r/t^/') (15) 

where P(r, t) is the probability that starting as nearest neighbors, the two probes are at a 
distance r apart at time t. In the range 1 < b < 2 this would imply that the average distance 
between the two probes grows as 

(r(t)) ~ (16) 

Since b is an increasing function of the Ising interaction e, this would predict a slower growth 
law of (r) with t, as e increases. 

The cumulative distribution function P(r, t) is defined as the probability that starting 
from a nearest neighbor position, the separation between the two probes at time t is larger 
than r. From Eq. [15] it follows that 

P(r,t) ~t(^-')/^r(r/t^/") (17) 

which means that P{r,t)t^''~^^^^ , plotted against r/t^/^ should show a scaling collapse for 
various values of t. 

In [12I, [1^ the time evolution of the average distance between the two probes was moni- 
tored numerically. Starting from a randomly disordered configuration, with the restriction 
that the two probes are placed on nearest neighbor sites, the system was evolved for a time 
tequii in an attempt to let it reach an equilibrium state. The time evolution during the equili- 
bration process followed the exchange rules shown in Eq. [T]with the important modification 
that the two probes were constrained to remain nearest neighbors i.e. they hopped together 
as if glued together. At the end of this equilibration, the medium is assumed to be locally in 
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steady state, in the vicinity of the probes, up to a distance of the order tg^^j;. At this point, 
defined as t = 0, the restriction for the relative position of the probes was released and 
the distance between them monitored. The distance between the probes was then assumed 
to follow the scaling form in Eq. [15] for t ^ tegmi when the two probes move within an 
equilibrated region. In this time regime, it was numerically verified that the growth of {r{t)) 
is consistent with Eq. [16] (3, Q- 

Note that the scaling form in Eq. [15] is expected to be valid in steady state. Therefore, 
to verify this scaling form, we followed the following procedure. Allow the system to reach 
steady state by evolving it without any restriction on the relative separation of the two 
probes. Then wait till the probes come to a nearest neighbor position with respect to each 
other and define t = at this point. Our data shows that (r(t)) follows Eq. [16] only for an 
initial time-regime, after which the growth exponent changes to ~ 1/3 which is close to the 
value of the growth exponent at e = 0. We present our data in fi^ 



We have also measured {r{t)) following the procedure of [12, [13(1 . We investigated the 
effect of different values of tequu and found the same behavior as described in the previous 
paragraph. Moreover, fig [1] shows that the curves for this partially equilibrated initial 
condition, coincide with that of the steady state initial condition (explained in the previous 
paragraph), for large time. 

We conclude that in steady state, {r(t)) does not follow Eq. [16] with an e-dependent b 
all the way, but shows a crossover at large time to the behavior t^^'^, which is the behavior 
obtained for e = 0. 



In [13| it was also reported that the cumulative distribution function P(r, t) shows a 
scaling form as in Eq. [TT] Starting from an initial configuration with the two probes next 
to each other (as discussed above), P{r,t) was numerically measured for a range of values 
of t and it was concluded that within that range, P(r, t)t^(''~^)/^ shows a scaling collapse for 
different values of t, as plotted against r/t^^"^. However, our numerical results indicate that 
this scaling collapse fails for larger t values (see fig[5lA-). Instead, an e-independent scaling 
form, more specifically, the scaling form expected for e = 0, seems to hold. We show this by 
plotting P(r, t)t^/^ against r/t^/^ and fig[5]B shows the scaling collapse for larger t values. 

In our simulation, we could not go to very large times as the finite size effects would 
become strong. In fig [5] we have presented our data for the largest system size {L = 4096) 
we could access. However, the crossover time is much smaller and no finite size effects are 
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FIG. 4: Average distance {r{t)) between the probe pair as a function of time. {r{t)) shows two 
different power law growths as time changes. The reference lines show that the growth exponent is 
(2 — h)/z at short times and changes to 1/3 at large times. The curves for partially equilibrated 
initial conditions (using the method of Rakos et al.) with different values of t^quU coincide for 
small t. We have also measured {r{t)) starting from steady state initial condition. The partially 
equilibrated data and steady state data coincide for large t. We have used e = 0.4 and L = 1000. 
Inset shows the steady state data for e = 0.3, 0.5 with L = 4096 and the reference line with exponent 
1/3. 



observed for this time-range. 

Our studies therefore show that the large time dynamics of the two probes is not affected 
by the presence of an interaction in the medium as reported in 



12 



13|, rather it resembles 



the case of non-interating medium. As discussed in Section 1, one possible rationalization 
is that the Ising measure of a KLS model (without probes) induces a finite e-dependent 
correlation length in the medium. In the presence of probes, the value of this correlation 
length may change, but it is expected to be finite still. As long as the displacement of 
the probes is less than this correlation length, the effect of varying e may be felt. But 
asymptotically, when the typical probe separation has exceeded the Ising correlation length. 
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FIG. 5: Color online: Panel (A) shows the lack of scaling collpse for the cumulative distribution 
function P{x,t) for larger values of t. As seen from the label of y axis, when P(x,t) is rescaled 
by an e-dependent prefactor, the scaling collapse works for smaller t, but fails for larger t. This is 
in contradiction with the scaling form described in Eq. Instead, P{x,t) is seen to follow a 
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it is plausible that they behave as if in a medium with no interactions, i.e. e = 22|. 



B. Macroscopic Number of Probes 

We now take up the study of a system with a macroscopic number of probes. We find that 
the dynamics of the STPs is governed by a diverging time-scale r, as in the non-interacting 
case e = 0. For t ^ r, an STP senses the fluctuations solely due to the KLS chain. But 
a KLS chain is known to have an Ising measure which means that if e is not too close to 
unity, only short-ranged correlations are present in the medium. Let tq be the time required 
for a probe particle to move a distance of order Rising- Then for tq <^ t <^ r, the dynamics 
of the probes in a KLS chain should be similar to those in an ASEP (where no correlation 

n 

is present in the medium) i.e. as that of the second class particles discussed in [6[. The 
dependence of the crossover time r on the probe density is discussed below. 

Let Ti be the separation between the i-th and {i + l)-th probe and Rm be the distance 
between the first and the (m + l)-th probe, i.e. Rm = Yl^i ^i- Let follow the distribution 
P(rj) ~ r^^. Assuming independence, the quantity Rm which is the sum of m such random 
variables should follow a Levy distribution with a norming constant ~ m^/'^^~^\ so long as 
Rm is less than the correlation length ^. In other words, the length Rm of a segment which 
contains m probes scales as m}/^^~^\ This is valid up to Rm ~ ^ but fails as Rm increases 
beyond that. ^ is the same correlation length that appears in Eq. [3] for the non-interacting 
case. Let m> be the number of STPs in a segment of length ^. Then m> ~ Hence in a 

system of length L, the total number of probes Nq can be written as A^^o = {L/^) which 
implies that the correlation length ^ ~ Po ^^"^ hence r ~ ~ Po^°^'^^"^\ where Zq is 

the dynamical critical exponent of the system. 

We have monitored the dynamical correlation functions Co(t), -Bo(^) and A(t), as defined 
in Eq. [5l [HI and [lOl respectively. Our numerical simulations indicate that these quantities 
follow the same scaling form as in the non- interacting case e = 6|. Moreover they continue 
to show crossover at a time-scale r ~ pg^, very similar to the e = case. In fig [6] we show 
the scaling collapse for Co(t) and 8^(1). We present our data for A(t) in fig [3 

In the case of two probes, one might expect A(t) would show the same scaling behavior as 
the second moment of the distribution P{r,t) in Eq. [151 ■'•S- A(t) should grow with time as 
our numerical simulations show that irrespective of the value of e, A(t) always 
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FIG. 6: Scaling collapse for Co{t) for e = 0.2 and po = 0.06,0.08,0.1,0.12. Inset shows scaling 
collapse for B{t) with e = 0.2 and po = 0.08, 0.1, 0.12, 0.15. We have used L = 16384. 
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FIG. 7: Scaling collapse for A{t) for finite e values. We have used po = 0.08,0.1,0.12,0.15 and 
L = 16384. The inset shows the linear growth of A{t) for L = 16384 in presence of two probes. 
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grows linearly with time (as with e = 0). We have shown our results for e = 0.5 in the inset 
in fig [71 

Note that the above scaling analysis and our numerical simulation presented in fig [6] and 



[7] point towards zo/(2 — A) = 3. If A = 6(e) as reported in [12], then for larger values of e 
this would lead to zq smaller than unity. For example, for e = 0.5, we have verified that the 
above scaling form remains valid (see fig[70, which would imply zq ~ 0.54 if A = 6(e). 

The other (simpler) alternative is that Zq = z = 3/2 and A = 3/2asine = case. This 
scenario would explain the observed po dependence of crossover time r. In the case of two 
probes, the above value of A is consistent with the large time growth exponent of the average 
separation (r(t)) between the probe pair (shown in fig[lj) and also with the linear growth of 
A(t) shown in the inset of fig [71 

V. CONCLUSION 

In this paper, we have studied the dynamics of shock-tracking probe particles in a one- 
dimensional KLS model of driven particles with nearest neighbor Ising interaction e. In 
particular, we have examined our results in the light of two different theoretical scenarios. 
The first scenario is based on an approximate mapping of the problem to a zero-range process, 
and leads to the conclusion that critical exponents characterizing power-law decays depend 
continuously on the strength of interaction, e. The second scenario is based on the premise 
that since the correlations induced by Ising interactions are short-ranged, asymptotic scaling 
properties which involve large distances and large times should be independent of e, and the 
same as at e = 0. The results of our numerical studies on the dynamical properties of the 
probe particles lend support to the second scenario. 

We find that in presence of only two probe particles in the system, starting from a steady 
state configuration where the two probes were nearest neighbors, the average distance {r{t)) 
between them shows a crossover in time. For an initial time regime (r (t)) the growth is 



consistent with a power law with an e-dependent exponent [2 — b{e)]/z |12l . [13| . However, 
for large enough time, the growth occurs with an exponent ~ 1/3, the value expected for a 
noninteracting medium, consistent with the second scenario discussed above. In addition, 
our study of the cumulative distribution of the probe-separation shows that for lar ge t ime, 
the distribution function does not follow an e-dependent scaling form as claimed in [l3| but 
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can be described by a form expected for e = 0, which again supports the second scenario 
mentioned above. 

For a small but finite density of the probes, the dynamical correlation functions show 
a similar scaling form as for e = [6]. These scaling forms involve a crossover time-scale 
r that diverges for small po as Pq^, as found for the non-interacting case We have 
seen that an e-dependent exponent 6(e) would lead to an e-dependent dynamical exponent 
Zq which may even become less than unity for larger values of e. The other option, an 
e-independent dynamical exponent zq = 3/2, is consistent with the second scenario outlined 
above, according to which, turning on a short-ranged Ising interaction in the medium does 
not change the large time and large distance properties of the system. 

It is not yet completely clear why the KLS-ZRP mapping does not seem to yield results 
which agree with the numerical results. One possible reason is that the ZRP mapping does 
not take into account the non-Markovian movement of the probes. This lack of agreement 
also opens up the question of the nature of the complete phase diagram for the problem 
under study, including negative values of e? 
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